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ABSTRACT 

Gamma ray burst (GRB) time series is a non-stationary time series with all its statis- 
tical properties varying with time. Considering that each GRB is a different manifes- 
tation of the same stochastic process we studied the time dependent as well as time 
averaged probability density function (pdf) characterizing the underlying stochastic 
process. The pdfs are fitted with Gaussian distribution function and it has been argued 
that the Gaussian pdfs possibly indicate the turbulence origin of GRB. The spectral 
and temporal evolution of GRBs are also studied through the evolution of spectral 
forms, color-color diagrams and hysteresis loops. The results do not contradict the 
turbulence interpretation of GRB. 
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1 INTRODUCTION 

Even after the three decades of extensive observations 
of gamma-ray burst (GRB) prompt emission there is 
no conclusive idea about the energy dissipation mech- 
anism and the radiative proces s responsib l e for high 
energy emission l|Lvutikovl l|2009l ): iGhisellinil (|2010l ')'). In 
the standa rd mod el of GRB, k nown as fireball model 
GoodmanI (Il986f): [Paczvnskil (1 19861 ): Shemi fc PiranI 
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■ .Sari fc PiranI (Il995l): iKobavashi et al.l (|l999h : 

iMeszaros fc ReesI (|20od '): iMeszaros et al.l l|2002l vr the bulk 



energy associated with the relativistic expansion of the 
fireball, is dissipated into the random energy of particles 
through internal shocks. Internal shocks are produced 
due t o the time depen d ent s pe ed of the relativi st ic out - 
flow (IKo bavashi et al.] (|l99it ): iRees fc Meszarod (|l994 ): 
ISari fc Piran (,199'^ )). But it is now established that the 
internal shock model suffers from different shortcomings, 
particularly, the low radiation efficiency, problem of baryon 
loadin g and the large dissipa t ion ra dius (|Naravan fc Kuman 
(|2009l ): lLazzati fc Begelmaiil (|2010l )'). 

As possible alterna ti ves of in ternal sh o ck m odel, 
iLvutikov fc BlandfordI (|2003h and iLvutikovl (|2006l ') pro- 
posed electromagnetic m odel of gamma-ray burst whereas 
iNaravan fc KumaJ (|2009l ') discussed relativistic turbulence 
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model for GRB. In electromagnetic model the GR B outflow 
i s con sidered as Poynting flux dominated wind (|Lvutikovl 
( 20061 )). The energy dissipation and particle acceleraton 
occur through magnetic dissipation process, mainly through 
magnetic r cconncction. In relativitic turbulence model by 
iNaravan fc Kumar (200i), it is considered that the fluid 
in GRB out flow is relativistically turbulent. The rapid 
variability in the GRB lightcurve arise due to the relativistic 
random velocity fluctuations of radiation emitting eddies 
formed in the turbulence. The model was successfully ap- 
plied to explain the X-ray-gamma-ray and optical emission 
during the prompt p hase of the brightes t Swif t gamma-ray 
burst GRB080319B |k umar fc NaravanI (|2009l )). 

A photospheric e m ission model is proposed by 
iLazzati fc BegelmanI (|2010l ) very recently. Generally 
the photosphere, associated with the flreball, is dominated 
by thermal radiation. It is assumed that the electrons, 
energised by some shock or due to magnetic reconnection, 
Comptonize the thermal photons in an optically thick 
medium. The multiple Comptonization leads to a power- 
law photon distribution in the keV-MeV region whereas a 
thermal peak remains in the keV energy range. The main 
advantage of of this model is that it does not require any 
efficient dissipation mechanism and the radiative efficiency 
of this model is also very high. But there is no strong 
observational evidence of any thermal component in GRB 
spectra. 

As far as the spectral studies are concerned, the majority of 
observed time averaged prompt emission spectra of GRBs 
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are fitted with Band function (|Band et alj (|l993l )) wiiich, 
indeed, is no more than a mathematical function having no 
physical interpretation yet. It is generally believed that the 
high energy emissi on in the prom pt phase is due to the syn- 
chrotron process (|Tavanil (|l996al lbh'). But the detail study 
of BATSE GRBs shows that for many GRBs the low en- 
ergy spectral slopes a re not consis t ent w ith the theory of 
synchrotron emission jPreece et alJ (|2000") ) . 
Therefore to have deeper understanding on the physical 
processes occuring in gamma-ray bursts it would be beni- 
ficial to study the time evolution of GRB spectra. It is 
also important to look into the GRB time series from a 
differ ent perspective and charact e rise th e underlying pro- 
cess. iLlovd-Ronning fc Petrosiaiil (|2002h : iGhirlanda et al.l 
l|2002h : iMedvedevI l|200^ and several other authors stud- 
ied the time evolution of GRB sp ectra for BATSE GRBs. 
iLlovd-Ronning fc PetrosianI l|2002l ) studied a 128 ms spec- 
tra for a large sample of GRBs observed by BATSE. The 
spectra were fitted with synchrotron emission model where 
the spectra were calculated in three emission regime, such 
as, for isotropic pitch angle emission of electrons, for small 
angle pitch angle distribution of electrons and for self- 
absorbed systems. Eventhough the spectral fitting were done 
and the time evolution of spectral index and normaliza- 
tion were studied, but the physical description was not 
self-co nsistently developed wit hin the framework of fireball 
model. IGhirlanda et al.l (|2002l ) also studied BATSE GRBs 
and fitted the time resolved the spectra with Band function, 
broken power-law, thermal Comptonization and synchrotron 
shock model. It was concluded that the behaviour of the time 
resolved spectra is not consistent with the time integrated 
spectrum. The low exergy spectral indices also violates the 
synchrotron spectral limit. 

In this paper we studied the probability density function 
(pdf) and spectral evolution of four bright GRBs detected 
by Swift and tried to find out if there is a ny sy stematics 
among those GRBs. In iBeloborodov et al.l (|200(]| '). the au- 
thors mentioned that the individual burst is a random re- 
alization of the same standard stochastic process. In any 
stochastic process there is an underlying pdf which charac- 
terises the process. Therefore it is important to find out the 
pdf from GRB time series to identify if it has any general 
characteristic. Since GRBs are highly non-stationary pro- 
cess, we calculated time-dependent as well as time averaged 
pdf for selected GRBs. We also studied the spectral evolu- 
tion, colour-colour and hardness-count rate correlations for 
those GRBs. 

In Section 2 we give a brief description of the selected GRBs, 
in Section 3 we describe the analysis procedure of Swift-'BAT 
data. Results are described and discussed in Sections 4 and 
we conclude the paper in Section 5. 
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Figure 1. First panel (from top) : Lightcurve of GRB080319b; 
second panel : probability density function - time averaged 
(black), time dependent (colours of different distributions cor- 
respond to the portions of the lightcurve with same colour) ; third 
panel : variation of spectral index with time; fourth panel : colour- 
colour diagram for the active portion of the lightcurve (same 
colour convention maintained). 



2 SELECTED GRBS 

We have analysed the Swift-BAT data for GRB050525, 
GRB061121, GRB080319b and GRB080411. GRBs are se- 
lected based on their fiuence. All four GRBs have fiuence 
> 10~^ erg cm~^. For the sake of completeness we give a 
brief description each GRB below. 

• GRB050525 : This GRB was observed with Swift- 
BAT in the energy range 15-350 keV. It ia at a redshift of 



0.606. The observed fiuence in the Swift-BAT energy range 
is 1.86 X 10~^ ergs cm^. The time averaged spectrum in 
the 15-350 keV range was fitted with a power-law (~ 
with spectral index 0. 83. The estimated i sotropic energy is 
Eiso = 1.94 X lO'^^ erg (|Band et all (|2005l ): lMarkwardt et al] 
(2005)). 

• GRB061121 : Swift-BAT detected this burst along- 
with a prescursor. Here we concentrate on the main event 
which has Tgo = 18.2 ± 1.1 s. The Swift/BAT spectrum 
was fitted with a power-law with spectral index 1.375. The 
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Figure 2. Same as Figurc[l] but for GRB050525. 
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Figure 3. Same as Figure[ll but for GRB061121. 



measured fluance is 2.792 x 10^^ erg cm~^. The measured 
redshift for the burst is 1.314 which le ads to a measured 
isotropic energy Eiso « 6.7 x 10^^ erg. iPage et al.l (|2007l ') 
studied the broadband spectral evoluton of the prompt and 
afterglow emission of the burst. 

• GRB080411 : GRB080411 consists of two strong 
peaks which are widely separated in time. The time aver- 
aged spectrum in the 5wift/BAT energy range was fitted 
with a power-law with spectral index 1.678. The observed 
fluence is 4.3 x 10~^ erg cm~^. The measured redshift of the 
burst is 1.03 which gives an isotr opic energy of the burst 
~ 4.04^^ erg (Mar shall et"ai (|2008i )'). 

• GRB080319b : This gamma-ray burst is one of the 
brighest GRBs with a very strong optical flash during its 
prompt emission. This GRB was simultaneously observed 
with Swift-BAT and if onus- Wind in the energy range 15- 



350 keV and 18-1160 keV respectively. The observed flu- 
ence in the 20keV-7MeV band is 5.7 ± 10"* erg cm"^. 
The measured redshift is 0.937 and isotropic energy is es- 
timated to be 1.3 X 10^* erg. The time averaged spectrum 
measured by ifonus-Wind was fitted by a broken power- 
law with break energy at Etr ~ 650 keV. The spectral 
index below 650 keV is 0.18 ± 0.01 and above 650 keV is 
2.87 ±0.44. This GRB was also followed by other telescopes 
in different wavelengths during the prompt and afterglow 
emission. During the prompt emission it was observed by 
the rORTORA telesc ope in the optical region and it is dis- 
cussed in detail bv iRacusin et al] (|2008[ l. This GRB i s stud - 
ied quite extensive l y by s everal authors (|Zou et al.l (|2009t ): 
iKumar fc NaravanI ()2009( )1. The most significant feature of 
this particular GRB is that the time averaged spectrum can 
not be explained by the synchrotron-self-Compton (SSC) 
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Figure 4. Same as Figurc[T] but for GRB080411. 



process in the framework of internal shock model. Rather 
a physically reasonable set of paramete rs can be obtained 
in the framework of turbulence model (iKumar fc NaravanI 
ll2009l) 1. 



3 DATA ANALYSIS 

We analysed Swift /BAT event data for four Swift detected 
gamma -ray bursts GRB050525, GRB061121, GRB080411 
and GRB080319b. The standard BAT software and the cal- 
ibration database (CALDB : 20090130) have been used to 
obtain the lightcurves and the spectra for all four GRBs. 
We extracted lightcurves with 10 ms time resolution in the 
energy band 15-150 keV for all GRBs. We also obtained 
lightcurves in the energy bands Ei : 15-25 keV, E2 : 25-50 



Table 1. GRB080319B: The fitted parameters, their uncer- 
tainities and the goodnoss-of-fit parameter for the Gaussian dis- 
tribution(s) corresponding to each pdf. The colours used to plot 
the pdfs are also given in column 1. For the first time segment, 
as the dof is zero, the fit is a representative one. 



Time 




Gaussian Parameters 






(sec) 


M 


A, 


Mi 




-10 4 


1 


429.7 


-0.07 


0.11 




(red) 












-4 — 10 




93.8 ±4.4 


0.66 ±0.03 


0.42 ±0.03 


23.16/22 


(green) 


2 


39.6 ± 3.0 


1.89 ±0.06 


0.40 ±0.04 


10 — 30 


1 


170.1 ± 4.6 


2.64 ±0.01 


0.44 ±0.01 


17.13/20 


(blue) 












30 — 52 




26.5 ± 8.6 


1.02 ±0.17 


0.23 ±0.09 




(magenta) 


3 


59.5 ±25.8 
184.8 ± 11.6 


1.57 ±0.08 
2.16 ±0.05 


0.21 ±0.1 
0.36 ±0.03 


27.46/18 


52 — 60 




235.9 ±44.9 


0.23 ± 0.004 


0.11 ±0.01 


0.85/3 


(cyan) 


2 


24.1 ±24.0 


0.43 ± 0.34 


0.21 ±0.13 


Table 2. GRB050525: Same as 


Table [T] For the last time seg- 




ment, as the dof is 


zero, the fit is a representative 


one. 




Time 




Gaussian Parameters 




■ xV'^ 


(sec) 


AI 


A, 






-3 — 


1 


96.5 ±9.5 


0.043 ±0.027 


0.086 ±0.014 


0.7/1 


(red) 












— 4 




20.6 ±2.2 


1.23 ±0.03 


0.32 ±0.03 




(green) 


3 


6.3 ± 1.4 
5.9 ± 1.0 


2.37 ±0.08 
3.82 ± 0.09 


0.28 ±0.08 
0.54 ±0.09 


40.8/37 


4—12 




34.0 ±3.5 


0.26 ±0.02 


0.19 ±0.03 




(blue) 


3 


29.9 ± 1.8 
17.9 ± 1.2 


1.24 ±0.04 
3.59 ± 0.03 


0.48 ± 0.04 
0.57 ±0.03 


28.5/41 


12 — 15 


1 


113.9 


0.001 


0.02 




(magenta) 













keV and E3 : 50-100 keV with time resolution of 200 msec 
and 1 sec for all GRBs. 

The spectra are extracted for every 5 sec from -10 sec with 
respect to the trigger time for GRB061121, GRB080411 and 
GRB080319b. In case of GRB050525 spectra are extracted 
for every 2.5 sec. 



4 RESULTS AND DISCUSSION 
4.1 Probability density function 

A stochastic process {x(t)}, also known as time series, is an 
ensemble of real or complex valued functions. Each stochas- 
tic process is characterized by a probability density function 
{pdf) depending on the underlying physic al process respon- 
sible for the generation of the time series ijBendat fc Piersoll 
(2010)). Therefore it is always interesting to develop an un- 
derstanding of the probability density function for a given 
time series. 

A time series can be stationary or non-stationary depend- 
ing on whether the statistical properties of the time series. 
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such as mean, variance, are independent of time or not. The 
gamma-ray burst time series is a real non-stationary time 
series with statistical properties of the time series vary with 
time. Here we calculate the time evolution of the pdf and 
the time averaged pdf for all four GRBs of our sample. To 
study the time evolution of the pdf we divide the GRB time 
series into different time segments and for each segment we 
calculate the pdf. If r be the time length of one segment 
starting at time t then the pdf is defined by 



Table 3. GRB061121: Same as Table[T] 



N 



(1) 



where AA'^i, the number of trials that the random variable 
falls between x and x -\- Ai, A'^ is the total number of tri- 
als. The error in p{x) at x is given by ^/AN^. For GRB time 
series x represents the count rate. The length of the time seg- 
ment r at time t is chosen through visual inspection of the 
GRB time structure. For the time-averaged pdf, p{x) is com- 
puted for the complete lightcurve. It is to be noted that in 
this work we have not used the normalization for p{x). The 
time-averaged pdf, pdfs for different time segments for each 
GRB, alongwith their respective lightcurves, are shown in 
Figures [T}{3] (/irst and second panel). The different portions 
of the lightcurve are shown in different colours and their 
corresponding probability density functions are also given 
with the same colour. The time resolution of the lightcurve 
used for the calculation is 10 ms for all GRBs. From the pdfs 
plotted for different time segment for each GRB it is evident 
that the pdfs corresponding to the rising and the falling part 
of GRBs peak around zero and those are not affected by the 
other portions of the GRB. As the GRB evolves, the pdf also 
evolves. In almost all cases, the pdf for each time segment 
is fitted with single or multiple Gaussian distributions given 
by 
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A - 








—3 — 10 


1 


408.1 ± 14.4 




0.14 ± 0.01 


1.6/2 


(red) 












10 — 60 


1 


2374 ± 144 


-0.01 ± 0.01 


0.1 ±0.004 


0.9/2 


(green) 












60 — 70 


2 


115.4 ±5.5 


0.52 ± 0.01 


0.26 ±0.01 


10.3/12 


(blue) 




49.5 ±4.2 


1.25 ±0.02 


n 1 s -1- n 09 
u. io m u.Ui 








54.8 ± 13.5 


0.35 ±0.02 


0.16 ±0.03 




70 — 80 




44.6 ± 2.8 


0.97 ±0.07 


0.39 ±0.14 


1 O c; /oo 
iZ.O/ 


(magenta) 


4 


37.2 ±6.1 
12.4 ±1.4 


1.85 ±0.05 
2.87 ±0.04 


0.24 ±0.04 
0.32 ± 0.05 


80 — 95 


1 


581.0 ±6.9 


0.09 ± 0.002 


0.10 ±0.002 


0.3/2 


(cyan) 












Table 4. GRB080411: Same as 


Table [1] 






Time 




Gaussian Parameters 






(sec) 
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A, 


Mi 




-10 — 


1 


359.2 ±26.8 


-0.02 ± 0.02 


0.14 ±0.01 


0.4/2 


(red) 












10 — 16 


1 


293.9 ± 18.7 


0.43 ±0.02 


0.20 ±0.01 


17.0/5 


(green) 












22 — 39 


1 


396.4 ±19.1 


0.08 ±0.02 


0.21 ±0.01 


8.3/5 


(magenta) 












48 — 70 


1 


442.3 ±23.5 


0.27 ±0.01 


0.20 ±0.01 


12.8/4 


(yellow) 













F{x) = ^ Gxp 



Mi 



(2) 



where M ^ 1. For all GRBs, the parameters of the fitted 
Gaussian function(s), their uncertainities and the goodness- 
of-fit parameter are listed in Tables [T}2] For the time seg- 
ments before the GRB trigger and after the GRB episode, 
the lightcurves are dominated by Poisson noise. Therefore 
the number of points obtained in the pdf for these time seg- 
ments are small. This results in very poor fit statistics. The 
variation of the parameters of the Gaussian distributions 
with the evolution of the burst reflects the non-stationarity 
of the underlying random process. We also attempted to 
fit the pdfs with log-normal functions but the fitting with 
Gaussian function was found to be statistically better than 
the log-normal distribution. There is no definite correlation 
between the peak position and width of the Gaussians. It is 
very clear that the time-averaged probability density func- 
tion profile (represented in black) can be fitted with multiple 
Gaussian functions. The possible physical interpretation of 
the Gaussian distribution is discussed below. It is argued 
that the Gaussian distribution possibly indicates the pres- 
ence of turbulence in the GRB outfiow originated in a rela- 
tivistically expanding fireball. Possible consequences of the 
turbulence in the spectral evolution of each of the selected 
GRBs are also discussed in detail. 



4.2 pdf and Turbulence 

It is shown by iKumar fc NaravanI (|2009l ) and IZou et ahl 
(1200911 that the optical and gamma-ray spectra in the 
prompt emission of GRB080319B can not be explained 
in the framework of internal shock model. Instead, 
iNaravan fc Kumaij (|2009l ) proposed the relativistic turbu- 
lence model to explain the prompt emission spectra of 
GRB080319B. In the framework of this model the highly 
variable gamma-ray lightcurve can be produced due to the 
emission from randomly moving eddies in the turbulence. 
The random relativistic variations of velocities of eddies 
cause the large amplitude variations in the individual pulses 
in gamma-ray burst lightcurve. In case of GRB080319B the 
optical emission is produced by synchrotron process in the 
inter-eddy medium whereas the gamma-rays are produced 
by inverse Compton process by energetic electrons in the ed- 
dies in the turbulence. To understand the probability density 
function obtained for the X-ray lightcurve of GRB080319B 
we argue as follows. 

If the large amplitude variations are due to the random ve- 
locity fiuctuations in turbulent eddies then the statistical 
properties of the velocity fluctuations are expected to be 
translated into th e random fiuctuati ons of the lightcurves . 
In recent works bv lMouri et"al] (|2002l ') and lChu et al] (|l996l ). 
it has been observationally shown that the probability den- 
sity function of turbulent velocity fluctuations is Gaussian 
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for fully developed turbulence. The probability density func- 
tion is sub-Gaussian when turbulence is developing and it 
is hyper-Gaussian during the decay of turbulence. We as- 
sume that these properties are also satisfied by the random 
velocity fluctuations in relativistic turbulence. In such a con- 
dition the intensity fluctuations in GRB lightcurve (which 
are supposed to be due to the velocity fluctuations of turbu- 
lent eddies) should also follow a Gaussian probability density 
function. In Figures [T}[4] (second panel) the pdfs are shown 
for GRB080319b, GRB050525, GRB061121 and GRB080411 
respectively. The pdfs for different time segments are shown 
in different colour and the time-averaged pdf is represented 
in black. The colour sequence is maintained same in all 
GRBs. It is to be noticed that the pdfs corresponding to the 
rising and falling edges of GRBs are completely separated 
from other pdfs corresponding to different time segments. 
In case of GRB080319b the pdf corresponding to a single 
segment is fitted with a single Gaussian function whereas 
for GRB050525 and GRB061121 pdfs corresponding to a 
single segment are fitted with multiple Gaussian functions. 
Multiple Gaussian distributions corresponding to a single 
lightcurve segment possibly indicate more than one region 
of turbulence in the GRB outflow. If we compare the rising 
and falling portions of GRBs with the development and de- 
cay of turbulence respectively then the corresponding peaks 
could be sub- and hyper-Gaussian distributions respectively. 
But such conclusion is very difficult to make with the present 
data. 

At this point it is important to mention that in presence of 
intermittency, t he turbulence stati s tics deviates from Gaus- 
sian behaviour (|Falcon et al.l (|2010l '): l\^ssilicosl (119951 )). Gen- 
erally when the amplitude of fluctuations are small, the tur- 
bulence statistics is Gaussian in case of fluid turbulence. 
In case of magnetohydrodynamic turbulence intermittency 
is generally present. The origin and nature of intermit- 
tency in fiuid as well as magnetohydrodynamic turbulence 
is not known. Therefore any departure of the peaks in ob- 
served probability density functions from Gaussian distri- 
bution does not rule out the present interpretation of pdf. 
Prticularly, if we consider the pdf of GRB080411 (Flgurel!]), 
we can see that the pdfs corresponding to the red, green, ma- 
genta and yellow portions of the lightcurve can be fitted with 
one or more Gaussian functions, but the pdfs corresponding 
to the peaks {blue and cyan) can not be fltted with Gaussian 
functions at all. It appears more like a power-law. Here we 
have not deliberately fitted them with any functional form. 
Having argued that the pdfs actually reflect the turbulence 
origin of GRB, now we concentrate on the spectral proper- 
ties of GRB to see if those properties can be explained in 
terms of turbulence. 



a non-thermal distribution. But GRB models in general do 
not discuss this issue in detail. In presence of turbulence 
which is presumably magnetohydrodynamic in GRB out- 
flow the particles can be accelerated by stochastic acceler- 
ation process or by magnetic reconnection process. In case 
of stochastic acceleration particles interact with the MHD 
waves developed in the plasma flow. This is a second order 
Fermi process and it is generally less efficient compared to 
the shock acceleration process. [Lazarian et al. (2003) dis- 
cussed that the stochastic acceleration process can be more 
efficient if momentum diffusion rate is higher than the pitch 
angle scattering rate and that is possible for a low density 
plasma in presence of high magnetic field. In such a con- 
dition the Alfven speed (in the units of c) exceeds unity 
and s uc h conditions ex i st in GRB outfiow (|Lazarian et al.l 
(|2003h ). iMao fc Wa"n3 l|201ll) discussed jitter radiation in 
small scal e random magnetic fiel d generated by turbulence. 
Following iHonda fc Hondal (|2005l ). they considered that par- 
ticles are accelerated in small scale random magnetic field. 
In such a condition, particle and photon spectra depend on 
the energy spectrum of turbulence. 

Otherwise also particle acceleration in presence of turbu- 
lence is possible through magnetic reconnection which is 
obvious in MHD turbulence. Turbulence possibly acceler - 
ates the reconnection process. iLazarian fc Vishniad (|l999D 
first put forward the turbulenc e reconnecton model and 
Ide Gouveia Dal Pino fc LazarianI (2003) discussed the par- 
ticle acceleration in turbulent reconnection. In general par- 
ticles can be accelerated almost instantly by magnetic re- 
connection process, but t his does not generate a power-law 
spectrum of particles. But 
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( 2003. ) and .Lazarian et al.l (|201ll ) showed that the particle 
acceleration in turbulent reconnection is a first order Fermi 
process and a particle's energy changes after each crossing 
of the reconnection zone. The average fractional gain in en- 
ergy is proportional to the reconnection speed. The magnetic 
fiux from two neighbouring regions of the fiuid approach 
each other with the reconnection spee d. It was shown by 
Ide Gouveia Dal Pino fc LazarianI (|2003l ) that the resultant 
particle spectrum is a power-law with spectral index -2.5. 
This is steeper than the spectrum (7~^) obtained with shock 
acceleration pr o cess, 7 being the particle Lorentz factor. 
iLarrabee etall (|2003l ) also studied the lepton acceleration 
by relativistic coUisionless magnetic reconnection. They ob- 
tained a distribution function for electrons which varies as 
~ 7~^, but the final spectrum (what the authors called ap- 
parent spectum) depends on the distribution function of the 
maximum Lorentz factor of electrons in each reconnection 
event. This makes the electron distribution ~ where 
g> 1. 



4.3 Turbulence, particle acceleraion and spectral 
evolution of GRBs 

4.3.1 Turbulence and particle acceleration 

The major challenge in high energy astrophysics is to under- 
stand the process of accelerating charged particles to rela- 
tivistic energies so that they can emit high energy radiation 
in a given physical condition. It is generally assumed that 
the first order Fermi acceleration across a shock (internal or 
external) accelerates particles to relativistic energies with 



4.3.2 Spectral properties 

We study the following spectral properties of GRBs of our 
sample to understand if turbulence can explain the observed 
results. 

• Spectral evolution : To study the spectral evolution 
of selected GRBs we extract the spectrum for every 5 sec for 
GRB080319B, GRB061121 and GRB080411 and for every 
2.5 sec for GRB050525. For the first three GRBs, the pho- 
ton spectra in the 15-150 keV energy range are fitted with 
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Figure 5. GRB080319B: Left panel : Hardness — count-rate dia- 
gram for energy bands Ei(15-25 keV) and ©2(25-50 keV), Right 
panel : Hardness - count-rate diagram for energy bands i?i (15—25 
keV) and E3 (50-100 keV). The colours used for different portions 
of the lightcurves (plotted as inset in each panel) are in direct 
correspondence with the colours used in hardness - count-rate 
diagram. 




Figure 6. GRB050525: same as Figure [3 



power-law (~ whereas for the last one the spectra are 

fitted with a power-law with exponetial cut-off. The spectra 
are fitted using minimization process within XSPEC. 
For GRB080319b, GRB061121 and GRB080411, the magni- 
tude of spectral index grossly increases as the burst evolves, 
indicating that the spectrum evolves from hard-to-soft with 
time. This is indeed observed for other bursts reported in 
the literature. The value of photon spectral index varies 
from ~0.6 to ~2.1 for GRB080319b, from ~1.3 to ~1.8 for 
GRB061121 and from to ~2.6 for GRB080411. It is 

to be noted that the evolution of spectral index is model 
dependent. Here the spectra are fitted with a single power- 
law which is statistically adequate due to the narrow en- 
ergy b andwidth of the detector. It was shown bv lBand et al.l 
1 19931 ) that the GRB spectra, over a broader energy range 
are well-fittted with a smoothly broken spectrum (Band 
spectrum). The observed spectral index variation reflects the 
evolution of underlying GRB spectrum which is intrinsically 
a smoothly broken spectrum with different spectral slopes 
on two sides of the break energy. If we consider the turbu- 
lent reconnection responsible for the acceleration of particles 
to relativistic energies and synchrotron process is the basic 
radiation mechanism then the value of photon spectral in- 
dex should be around 1.75 for instantaneous spectrum and 
2.25 for cooling dominated spectrum. But the value of ob- 
served photon spectrum indices fall below these limits. This 
implies that even if turbulent reconnection plays an impor- 
tant role in accelerating particles in GRB outflow, it is not 
the only mechanism working there. This requires further de- 
tailed study. 

In case of GRB050525, the spectral evolution pattern is not 
very clear because of its short duration. The spectra of this 
GRB could not be fltted well with power-law, instead a 
power-law with exponential cut-off fitted the spectra well. 



• Colour-colour diagram : The hardness ratio between 
the two energy bands, say Ei and E2, is defined as 

HR,, = g (3) 

where Ci and C2 are the count rates in the energy bands Ei 
and E2 respectively. The error in the estimation of hardness 
ratio is 

Similarly the hardness ratio and its error AHR13 can 

be defined for the energy bands Ei and E3. 
Colour-colour diagrams, shown in Figure [T}{4](/onri/i panel), 
show the correlations between the hardness HR12 and 
HR2i. These plots also have similar colour sequence as the 
lightcurves and are generated for the active part of the GRB 
lightcurvers. The time resolutions used for these plots are 1 
sec (GRB080319b, GRB080411) and 200 msec (GRB050525, 
GRB061121). It is generally found that the colour-colour 
diagrams show systematic correlations over the complete 
evolution of all bursts. The striking feature for the colour- 
colour diagrams for all GRBs is that they start with a very 
hard spectrum. This is possibly due to the almost instanta- 
neous acceleration of particles to high energies. This could 
be possible in magnetic current filaments generated in turbu- 
lence. T he estimated accelera tion time scale for such event is 
given bv lMao fc mn3 ()201ll ). It depends on magnetic field 
strength and turbulent length scale in the fiow. Later part of 
the evolution in colour-colour plane is primarily dominated 
by the cooling of particles. This makes the particle spectrum 
steeper and so the photon spectrum. It is important to note 
that in case of GRB080319b, the blue and magenta points 
(Figure [1] fourth panel) in colour-colour plane are clustered 
in a small region. These points correspond to the nearly 
flat portion of the lightcurve (represented in same colours). 
During this portion particle acceleration rate and cooling 
rate were comparable, thereby making the hardness almost 
constant. 

• Hardness— Count rate correlation : The hardness 
and count rate correlations have been studied for the ac- 
tive phase of each GRB for the energy bands i?i(15-25keV), 
£2(25-50keV) and £;3(50-100keV). For each burst, hard- 
ness vs countrate plots are shown in Figures [5] - |8l The left 
panel in each figure shows the hardness vs count rates cal- 
culated for energy bands Ei and E2 while the right panel 
shows the same for energy bands E\ and E^,. The repre- 
sentative colours in the lightcurve show the corresponding 
portions mapped into the hardness-count rate plane. The 
lightcurves used for this calculation had a time resolution of 
200 ms except for GRB080319B where the lightcurves had 
time resolution of 1 sec. All hardness-count rate loops for 
GRB080319B, GRB080411 and GRB050525 are in clockwise 
sense except for the burst GRB061121. For GRB061121 the 
loop structure is quite complex. It starts with an anticlock- 
wise loop for the first peak {red portion), then it traverses 
clockwise path for the second {blue) and it goes anticlock- 
wise for the third peak {pmk). The clockwise loops repre- 
sent softjagwhereas the anticlockwise loops indicate hard 
lag (iKirk et al.l l|l998l )). This is possibly due to the fact that 
during the fast rise of the spike, more particles are acceler- 
ated quickly to higher energies by some acceleration process. 
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Figure 7. GRB061121: same as Figureg] 




Figure 8. GRB080411: same as Figure[5l 



During the decay part of the spike the coohng of the particles 
dominates and generates soft lag in the spectral evolution. 
This feature is very well reflected in the hardness-count rate 
loops shown in Figures [5]-[8l The clockwise loops arise when 
the evolution of spectral slope is primarily governed by any 
cooling pro c ess w hich is faster for higher energy particles 
ijKirk at all l| 19981 )). For GRB061121 anti-clockwise loops, 
signifying hard lag, appear when the acceleration time scale 
is nearly equal to the cooling time sca le and part i cles a re 
gradually accelerated to higher energies l|Kirk et al.l (| 19981 )). 
These behaviours of spectral evolution of GRBs strongly de- 
pend on the particle acceleration process and the radiation 
emission mechanism occuring in the GRB outflow. With the 
present understanding of particle acceleration process in tur- 
bulence, it is not completely understood what kind of parti- 
cle acceleration and cooling process can really explain such 
properties along with the evolution of the exact spectral 
forms. 



5 CONCLUSION 

In this paper we studied the time-dependent probability 
density function obtained from the observed lightcurves of 
four bright Swift gamma-ray bursts. The probability den- 
sity functions are fitted with Gaussian distribution function. 
The parameters of the distributon function vary with time 
as the burst evolves. Considering the fact that the short 
time-scale fiux variations in GRBs are due to the random 
motion of radiation emitting turbulent eddies, as proposed 
bylNarayan & Kumar (2009), we argued that the Gaussian 
probability density function actually refiects the presence of 
turbulence in GRB outflow during the prompt emission. 
We also studied the spectral properties of GRBs, such as the 
time evolution of spectral index, colour-colour diagram and 
hysteresis loops to understand if they can be explained in 
the framework of turbulence. With the present understand- 
ing of particle acceleration in turbulence, it is found that 
the particle acceleration process in turbulent recconection 
can partially explain the observed spectral indices in 15-150 



keV energy range. The behaviour reflected in colour-colour 
diagrams are consistant with spectral evolution seen for all 
GRBs. But deeper understanding of colour- colour diagrams 
and hysteresis loops require detailed knowledge of different 
time scales, therefore a detailed study of radiative processes 
in GRBs is essential and is beyond the scope of this paper. 
Finally, we would like to rem ark that for steady sources like 
X-ray binaries (e.g. Cyg X-1, (Wu et al.' (2010'); Uttlcv ct ajj 
((2005 ))) and flaring systems like solar flares ( Zhang (2007,))) 
the probability density functions are found to follow log- 
normal distributions. For both types of systems turbulence 
seems to a major role. In the light of such observations the 
main finding of this work, i. e. the Gaussian probability den- 
sity functions for GRBs needs further study. 
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